Theory of proximity effect in normal metal/ d^2_y2-wsive superconductor interface in 
the presence of subdominant components of the pair potentials 
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Superconducting proximity effect in normal metal (N) / d^2_j,2-wave superconductor (D) junctions 
in the presence of attractive interelectron potentials which can induce subdominant s-wave pair 
potentials both in N and D sides, is studied based on the quasiclassical Green's function theory, 
where spatial dependencies of the pair potentials are determined self-consistently. In the N/D 
junctions with orientational angle with 6 = 0, the s-wave component is induced in the N side by 
the proximity effect only for high transparent case, where the induced s-wave components in both 
the N and D sides do not break the time reversal symmetry (TRS). For fully transparent case, the 
resulting local density of states has a very sharp zero-energy peak (ZEP), the origin of which is the 
sign change of the pair potentials felt by the quasiparticles between the s-wave component in the 
N side and d^2_y2-wave dominant component in the D side through Andreev reflections. On the 
other hands, for = n/4:, the subdominant s-wave component which breaks the TRS appears near 
the interface. Besides, for lower transparent cases, the subdominant imaginary s-wave component 
is also induced near the interface in the N side. The proximity induced s-wave component in the 
N side does not enhance the magnitude of the s-wave component of the pair potential which break 
the TRS in the D side. The resulting LDOS at the interface has the ZEP or its splitting depending 
on the transparency of the junction. 

PACS numbers: 74.50.-|-r, 74.20.Rp, 74.72.-h 



I. INTRODUCTION 



To determine the pairing symmetry in unconventional 
superconductors is an interesting problem to understand 
the pairing mechanism of superconductivity. Nowadays, 
it is widely accepted that the superconducting pair po- 
tential of high- Tc cuprates has a dx'^-y2-wa,Ye symmetry 
in the bulk stateiSiii^ In order to determine the pair- 
ing symmetry, several phase-sensitive probes have been 
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Among them tunneling spectroscopy via An- 
dreev bound states (ABS's)|^^i^i'^d°did^d3i^^d^i^^ help 
us to detect the internal phase in the pair potential. The 
formation of ABS's at the Fermi energy (zero-energy), 
which is originated from the interference effect in the ef- 
fective pair potential of dx2^y2--wave symmetry through 
the reflection at the surface/interface, plays an impor- 
tant role when the angle 6 between the lobe direction of 
the c?2,2_j,2-wave pair potential and the normal to the 
surface/interface is nonzero.^ In fact, the existence of 
ABS's, which manifests itself as a distinct conductance 
peak at zero-bias in the tunneling spectrum (ZBCP), has 
been actually observed. Up to now, the consistency be- 
tween thcoric j^^'^^'^^'^"'^^'^'^'^'^'^"^'^^'^^'^'''^^'^^'^"''^^'^^''^^''^^ 

and experimenta35,36i37,38,39,40,41,42,43,44,45,46,47,48,49 

been checked in detail. 

On the other hand, at the surface/interface of (i^2_j,2- 
wave superconductor, it is known that the time-reversal 
symmetry (TRS) may be broken The reduction 
of the amplitude of (i^2„j^2-wave pair potential near 
the surface/interface at low temperatures allows to in- 



duce the subdominant s-wave [dxy-'wave] component in 
the imaginary part of the pair potential, i.e., d-l-is- 
wavo"d^i^^i^°i53,54 [d-^id'-wave^^i^^ ] state. The am- 
plitude of the induced s-wave component has a maxi- 
mum at 9 — 7r/4, which blocks the motion of quasi- 
particles. Then, the energy levels of ABS's shift from 
zero and the resulting tunneling spectra has a ZBCP 
splitting even in zero magnetic field. The observation 
of the ZBCP splitting without magnetic fields was be- 
lieved to be one of the evidence for the broken time re- 
versal symmetry states (BTRSS). However, actual tun- 
neling experiments still remain to be controversial. In 
fact, some groupa39i4°i^3,44,45,46 j^^^^g reported the ZBCP 
splitting and they ascribed the origin of the ZBCP split- 
ting to the above BTRSS. At the same time, there are 
other experiments which do not show the ZBCP split- 
ting, and in these experiments, the ZBCP survives even 
at low temperaturei37i38,4i, 42,47, 48,49 Moreover, a criti- 
cal current measurement of grain boundary junctions 
in high-Tc cuprates concluded the absence of BTRSS 
at the interfaceiS It has been studied in detail that 
the ZBCP splitting in the tunneling spectra is sensitive 
to several factors: (i) transmission probability of the 
junctionSfS (ii) roughness at the interface and 
(iii) effect of impurity scattering in normal metals or 
superconductorsi ^3,25 Recent studies by Asano et al.^ 
revealed that the existence of the impurity scattering can 
induce the ZBCP splitting even without BTRSS. Taking 
account of these situations, one can not conclude that 
the observation of the ZBCP splitting by tunneling ex- 
periments directly indicate the presence of the BTRSS 



near the surface/interface. 

More recently, it is proposed that the induced s-wave 
component of the pair potential by proximity effect in 
normal metals (N) may enhance the magnitude of the 
subdominant s-wave component in dj.2_j^2-wave super- 
conductors (D), which forms the BTRSS based on the 
analysis of the tunneling experiments.®*' It is of great in- 
terest to study whether the induced s-wave component in 
the N side by the proximity effect has the influence on the 
subdominant component in the D side which forms the 
BTRSS. The proximity effect in the N/D junction with- 
out the BTRSS was theoretically studied by Ohashi.®-'^ It 
is shown that the amplitude of the subdominant is-wave 
is the largest at 9 — tt/A. For 9 = Q with fully transparent 
junctions, local density of states (LDOS) at the interface 
has a zero energy peak (ZEP), the origin of which is not 
ABS but the sign change of the pair potentials between 
N and D sides felt by Andreev reflection of quasiparti- 
cles. It is necessary to study the interplay between the 
proximity effect and the BTRSS by changing the orienta- 
tional angle 9. At the same time, it is also interesting to 
clarify the relationship between the different origin of two 
kinds of ZEP's. One is the ZEP which originates from 
the formation of ABS for 9 ^Q. The other is due to the 
proximity effect expected for fully transparent junctions 
for 61 = Oi^ 

In the present paper, we study the proximity effect in 
N/D junctions on the basis of quasiclassical Green's func- 
tion methods. We assume attractive interelectron poten- 
tials which induce subdominant s-wave components in 
the N side as well as the D. The spatial dependencies 
of the pair potentials both in the N and D sides are de- 
termined self-consistently. The LDOS at the interface of 
the N/D junctions are studied in detail by changing 9 and 
the transparency of the junction. For 9 — 7r/4, the mag- 
nitude of the s-wave component is induced even in the 
N side, while the subdominant s-wave component which 
breaks TRS exists near the interface on the D sidc.*^^ On 
the other hand, for 9 — Q, s-wave component is induced 
in the N side by the proximity effect with the increase 
of the transparency of the junction. The induced s-wave 
component in the N side and that in the D side do not 
break the TRS. It is revealed that the proximity induced 
s-wave component in the N side does not enhance the 
BTRSS in the D side. 

The organization of the present paper is as follows. In 
Sec. ^ a theoretical formulation to calculate the spatial 
dependence of the pair potential, and local density of 
states is presented. In Sec. IIIII results of the numerical 
calculations are discussed in detail. Finally, we summa- 
rize the paper in Sec. IIVI 



II. THEORETICAL MODEL 

In the present paper, we consider the N/D junction 
separated by an insulating interface at x = 0, where the 
normal metal is located at a; < and the (i„2_„2-wave 




FIG. 1: Schematic illustration of a N/D junction with a spa- 
tially varying pair potential. The misorientation angle 9 de- 
notes the angle between the crystal axis of ci^2_j^2-wave su- 
perconductor and the normal to the interface. 



superconductor extends elsewhere (Fig. ^ . For the sim- 
plicity, two dimensional system is assumed and the x axis 
is taken perpendicular to the interface. When quasiparti- 
cles are in the xy plane, a transmitted electron like quasi- 
particle and holelike quasiparticle in the D side feel dif- 
ferent effective pair potentials (</>+) and (</()_), with 
0+ — <f> and (p- ~ TT — (j). Here (j) is the azimuthal angle 
in the xy plane given hy {k^ + iky)/ \k\ — e"^. Besides, 
a cylindrical Fermi surface is assumed and the magni- 
tude of the Fermi momentum and the effective mass are 
chosen to be equal both in the N and D sides. 

In order to study the proximity effect in the N/D junc- 
tion, we determine the spatial variation of the pair po- 
tentials self-consistently. For this purpose, we make use 
of the quasi-classical Green's function procedure^ devel- 
oped by Nagai and co-workers In the following, we 
briefly summarize this scheme we employed. We intro- 
duce the generalized Eilenberger equationjS 
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(1) 
(2) 



where upx = vpcos9 and fi{i — 1,2,3) stand for the x 
component of the Fermi velocity and the Pauli matrices, 
respectively. Here ujm = T:T(2m + 1) (m: integer) is the 
Matsubara frequency. 

Now, considering a semi-infinite N/D junction geome- 
try, the pair potential in D [N] side tend to the bulk value 



zero 



A 



[A 



-oo)] at sufficiently large x. In 



semi-infinite limit, we can find the quasi-classical Green's 
function g^ai't'a.x) for l{= N,D) regions given bjiSMiSfe 



2AL(:^) 
TrKW: 



(3) 
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Following equations are satisfied for the N side, 

- J7^(0+,:E,O)i?DA^(O,:E), (4) 
A^{x) ^ \^ix,0)RBU^{(t)^,0,x), (5) 



X^{x,x)—i •^N/^^ _r.J^^^\ /r.N{^'\ I' 



and for the D side, 



= AD(x,O)i?N?7°(0+,O,a:), (6) 
A^{x) ^U^{(j)-,x,0)RNX^{0,x), (7) 

ig^ix) -g^{x)/g^{x') 



respectively. In the above, i?N and i?D are matrices with 
reflection probability R at the interface given bySi 



«Noci i(i_i^)^N(Q)^N(o) gl{Q)R-g^{Q) 
1^ i(i - i?)^D(o)eP(o) gD(o)i? _ 00(0) J • 



Here, we define g^ix) given by 



Un^^\(t)a,x) 



Un'' \4>a,x) 



(10) 



with 



(15) 



(XN(a;) FN(^))^^^(^N(+)(^^^Q^ «r+)(0+,O) Jr2 
X i?D?7!'(</'-,0,x), (16) 
(XD(X) yf(a;))f2 = (^,°(-)(0_,O) ^;,?(-)(0_,O))f2 



X i?N?/+(0+,O,x), (17) 



(0+,o) ; 

(18) 



Then we define ^a{x) given by 



(19) 



Next, in order to obtain the quantities (7^(0q,, a;, a;') in 
Eqs. ©-(Pi, we rewrite ia(x) asS^ 



T2, 
(11) 



i^(x) = [ ] ( ^^(^) ) (12) 



Vn''^'{(t>+,x) 



A^^{x)^{ V) !; 1 (^?(-) (13) 



X^ix) 
Y^{x) 



_,X) Vn'- \4'-,x} J T2, 
(14) 



Here, ^/^(x) and J^^(a;) in Eqs. (0 and (O, obey the 
following Riccati type equations: 



n\vF.\-g^gUx) =a [2Lu„,gUx) + A'{(t>a,x)gUxf 



A\<p^,xr] , 

+ A'(0„,x)] . 



?i|«F.|^^i(x) = - a [2w™.Fi(:r) - A'(,/)„, x)* J-^(x)2 



(20) 
(21) 



We can write the quasi-classical Green's function in a 
compact form— 



._^^^_/^l + e^(x)^^(-) 2i^N(^) \ 

5D . - ^^ _ i ( + 1 2i.FD(x) \ 

g++(<^+,x)- 2igD(^) -e°(x)^°(x)-l j' ^^^^ 

I 



with g_ {(j)^,x) — —5^ I (—(/>+, x)T. Initial conditions of these equations are as follows. 



A°(0„,ai)' 



5_»(-oo)=0, 50,„, = _^:^, ,24) 
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FIG. 2; (a) Spatial dependence of the pair potentials in a 
normal metal /s-wave superconductor (N/S) junction with 
T = 0.02Ts and Tn/Ts = 10"^. (b) The corresponding local 
density of states (LDOS) at the S side of the interface. 



J 12, 



(29) 



where luc is the cutoff energy and [5^^ {(pa , a;)] 12 means 
the 12 element of 5^„((/)a, a;). Here V^tji^cjia) is the ef- 
fective inter-electron potential of the Cooper pair in the 
I side. In our numerical calculations, new A'(0a,x) and 
Qaai^a^x) are obtained using Eqs. (1201-1(221) and (j^HJl. We 
reiterate this process until the convergence is sufficiently 
obtained. 

Based on the self-consistently determined pair poten- 
tials, the LDOS can be calculated as, 



1 r'^ 

Ni{E,x)^- d(t>ni{E,(t),x), 

•^J-7r/2 



(30) 



*E+i5 



(31) 



where Nq means the density of states (DOS) in normal 
states, and S is infinitesimal. In this paper, we choose the 
temperature T as T/Td = 0.02, where Td is the critical 
temperature of the bulk c?2,2_j,2-wave superconductor. 



III. RESULTS OF NUMERICAL 
CALCULATIONS 



with = \/uj^ + |A'((/)a, oo)p. Moreover, the bound- 
ary condition of the ^^(0) and 7-^(0) at the interface 
a; = are 



In this section, we show our results of numerical calcu- 
lations on the spatial dependence of the self-consistently 
determined pair potentials and the corresponding LDOS. 
The spatial variation of the pair potentials in the N/D 
junction is expressed as 
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(25) 
(26) 
(27) 
(28) 



The pair potentials for both N and D sides are given 



AN((/.,a;) = AN(x), 
/S°{<i),x) = Ad{x)cos2{(j>- 



(32) 

9) + A,{x), (33) 



where AN(a;), Ad{x), and As{x) denote the amplitude of 
s-wave component in the N side, d^2_y2-wave component 
in the D side, and subdominant s-wave component in the 
D side, respectively. The attractive potentials V''{(l),(j)') 
with / = N, D, are given by 



^^(0,0') 



2VdCos(2(/)- 26*) cos(20' 



29) + Vs., 



(34) 
(35) 
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FIG. 3: Spatial dependence of the pair potentials in N/D 
junctions with 6 = 7r/4 and T^/Td = 0.01. (a) ii = and 1, 
(b) R = 0.25, 0.5, and 0.75. 



and 
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Td ^ TO + 1/2 



27rfcBr 



(36) 
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Here, Tn denotes the transition temperature of s-wave 
pair potential in the N side. Tg denotes the transition 
temperature of s-wave pair potential in the D side with- 
out dx^-y2-w&ve attractive potential. 



First, in order to understand the role of the proximity 
effect clearly, let us check the case of normal metal/s- 
wave superconductor (N/S) junctions. FigureEIa) shows 
the obtained spatial dependence of the pair potentials in 
the junctions for extremely low {R — 1) and high {R = 0) 
transparency cases. The x-axis of Fig.[21is normalized by 

= ^"^^f/ttAso, which is the coherence length of s-wave 
superconductor. For i? = 1, the s-wave pair potential re- 
mains to be constant in S side, whereas the pair potential 
in the N side is zero. The spatial variation of the pair 
potentials is represented as a step function. The result- 
ing LDOS at the interface reproduces the bulk U-shaped 
DOS [see Fig. [Hb)]. On the other hands, for i? = 0, 
the pair potential An in the N side survives toward the 
inside. We can see that the LDOS at the interface in the 
presence of the proximity effect is different from the bulk 
DOS. 



A. N/D junction with 6 = 7r/4 

In this subsection, we focus on the N/D junction with 
(llO)-oriented interface {6 = it/ A). As the parameter of 
the subdominant s-wave component, we take Ts/Td — 
0.3 in the following. In Figs. |2Ia) and|2Ib), the spatial 
variations of the pair potentials is shown for T^/Td — 
0.01 for i? = 0, 1, and 0.25, 0.5, 0.75, respectively. Here, 
Re[A(;(a;)] and Im[AN,s(2;)] denote real part oi Ad{x) and 
imaginary part of AiM,s(a;), respectively. The x-axis of 
Fig. 13 is normalized by the (i^2_j,2-wave coherence length 
^rf — hv-p/nAdo, where A^o is dx2_y2-w&ve pair potential 
in bulk states. 

First, we present our results in the light of previ- 
ous theories liSiSS The reduction of Re[Ad(x)] originates 
from a depairing effect that the effective pair potentials 
A°(0-|.,O) and A° ((/>_, 0) have reversed contribution to 
the pairing interaction for certain range of </>. Re[Ac;(2;)] 
is suppressed at the interface in the D side. At the same 
time, the quasiparticle forms the ABS with zero-energy 
at the interfaceiSS When the magnitude of the reflection 
probability R approaches unity, the ABS becomes unsta- 
ble with the introduction of the subdominant s-wave at- 
tractive potential in the D side. And then Im[As(a;)] is 
induced in the vicinity of the interface in the D side in 
low transparent cases with large magnitude of Rm^ As re- 
gards the proximity effect in the N side, the pair potential 
An(x) is not induced for fully high and low transparent 
cases. However, the imaginary component of An (a;) is 
enhanced nearly in i? = 0.75 [see Fig. We can rec- 
ognize that the existence of Im[As(a;)] near the interface 
of D side can allow the enhancement of the Im[AN(a;)]. 
This fact is more recently found by Lofwander,^^ and our 
results are consistent with his work. 

In Figs. Ela) and E^b), we show the corresponding 
LDOS at the D side of the interface, where the same pa- 
rameters of Figs.|3Ia) andOIb) are used, respectively. For 
i? = 0, N-[){E, 0) is equal to A'o and has no E dependence. 
If we neglect the induced s-wave component in the N 
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FIG. 4: The LDOS for the N/D junction with Q = 7r/4 and 
Tk/Td = 0.01. (a) ii = 0, 0.25, and 0.5, (b) R = 0.75 and 1. 
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FIG. 5: Spatial dependencies of the pair potentials in N/D 
junctions with = and T^jTi = 0.01. (a) il = and 1, (b) 
R = 0.25, 0.5, and 0.75. 



the previous worksiii^^iS. 

B. N/D junction with 6 = 



side, the resulting LDOS has a zero-energy enhanced line 
shapes for small i?.^^ This is because that the magnitude 
of the subdominant s-wave component, i.e., Ini[As(a;)] 
is very small and the incident and reflected quasiparti- 
cles normal to the interface can feel opposite signs of 
the dj.2_j^2-wave pair potentials. However, as shown in 
Fig-EJb), not only Im[AN(a::)] but also Im[As(x)] remain 
to be non-zero due to the proximity effect. In this case, 
the resulting LDOS has complex line shapes. For i? = 1, 
due to the enhancement of the magnitude of lm[As(2;)], 
the LDOS has the ZEP splitting, which is consistent with 



In this subsection, we focus on the N/D junction with 
6* = and T^/T^ = 0.01. The spatial dependencies of the 
pair potentials Re[AN(a:)], Re[Ad(2:)], and Re[As(a:;)], are 
plotted in Figs.EIa) andEIb), for i? = 0, 1, and 0.25, 0.5, 
0.75, respectively. For the low transparent limit, R — 1^ 
[see Fig.ini^a)], the amplitude of the pair potential is con- 
stant. This situation is similar to the case of the N/S 
junctions with i? = 1 [see Fig. |2Ia)]. With the decrease 
of R, Re[Ac;(2;)] is suppressed near the interface, while 
Re[As(a;)] is slightly mixed in the vicinity of the inter- 
face. Since Im[As(a;)] = is satisfied, the TRS is not bro- 
ken near the (100) interface. On the other hand for the N 
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FIG. 6: (a) The LDOS at the interface by its value in the 
normal state for various R with ^ = and T-^/Td — 0.01. (a) 
The N side, and (b) the D side. 



side, we can readily see that the pair potential Re[AN(2;)] 
is induced. This indicates that the superconducting pair 
potential penetrates into the N side due to the proxim- 
ity effect. With increasing the magnitude of Tn, the 
amplitude of Re[AN(a;)] becomes larger. However, both 
Im[As(a;)] — and Im[AN(a;)] = are satisfied. This 
situation is significantly different from the corresponding 
case oi 6 — tt/A where the TRS is broken. 

Next, let us look at the corresponding LDOS at the 
interface. We show that the LDOS at the interface in the 
N and D sides as shown in Fig.El For R= 1, the LDOS 
on the D side has a V-shaped structure similar to the 
bulk d-wave density of statesiii*SiS The corresponding 
LDOS in the N side is constant. With the decrease of 
the magnitude of R, the magnitude of LDOS in the D 
side around zero energy is enhanced while that around 
±Ado is suppressed. For R = 0.25, the LDOS has a 
small dip like structure both in the N and D sides. The 
fine structure around zero energy is due to the induced 
pair potential An (a:) in the N side. The extreme case 
is i? = 0, where normal reflection is absent. The pair 
potential Re[AN(a;)] induces the ZEP in the LDOS [see 
Fig. |S] through Andreev reflectioniSi 

Now we concentrate on the width of this ZEP. We see 
that the relevance of the peak width and the infinites- 




FIG. 7: The LDOS at the interface near zero energy with 
R = for various S. 



imal number 5 introduced in Eq.|j25 in order to avoid 
divergence in the actual calculation where the inverse of 
S can be regarded as a life time of quasiparticles. As 
shown in Fig. |3 the width of the ZEP becomes narrow 
with the decrease of the magnitude of S, while the height 
of it increases monotonically. For S ^ limit, the ZEP 
is reduced to be expressed by the (5-function. 

In Fig. IHI we show the spatial dependence of LDOS in 
the N side {x < 0). For i? = 0, the zero energy states 
formed at the interface penetrate into the N side. The 
width and the height of the ZEP is reduced with the in- 
crease of the magnitude of x. For R = 0.1, the LDOS 
has a mini gap. The width of the gap has a spatial de- 
pendence since the induced pair potential Re[AN(2;)] de- 
pends on X. The width of the mini gap is reduced with 
the increase of the magnitude of x. 

At the end of this subsection, let us remark on the 
difference on the origin of ABS's between 9 = tt/A and 



9 = 0. As mentioned in Sec. lIII Al the formation of ABS's 
is originated from the opposite sign of the pair potentials 
between Ad(0+,O) and Ad{(l)-,0) felt by quasiparticles 
in the D side. However, in the case oi 9 = 0, there is 
no sign change between Ad{(t)+, 0) and A^ ((/)_, 0) for any 
(f). This manifestation of the ZEP in the LDOS is due 
to the another origin. It is due to the sign change of the 
pair potentials between Re[AN(a;)] and Re[As(a;)] or that 
between Re[AN(a;)] and Re[Ad(a;)] for |0| > tt/A through 
the Andreev reflection. Actually, by calculating the angle 
resolved LDOS nD(£', (/), x) in Ea. (|^ . n£,{E,(j),x) has 
the ZEP for cj) — tt/3 while it does not for cj) — tt/6 as 
shown in Fig. |51 




(b) R = 0.7 
0=0 




FIG. 8: The LDOS at a; < with 6 = 0. (a) 7? = and (b) 
R = 0.1 for various x. Here = vf/2hT is the coherence 
length in the N side with T = 0.02rd. 



C. N/D junction with 6* = 7r/6 

In this subsection, we concentrate on the case with 
intermediate angle 0, e.i?., 9 = 7r/6. The spatial depen- 
dencies of the pair potentials become very complicated 
as shown in Fig. E| where we choose T^/Td = 0.01 for 
various R. For i? = 0, the amplitude of Re[AN(a;)] is 
induced toward the inside of N. This situation is simi- 
lar to the case of 6* = by the proximity effect. With 
the increase of R, the amplitude of Re[Ad(a;)] is signifi- 
cantly suppressed at the interface, due to the destructive 
interference originating from the sign changing nature of 
da;2_j,2-wave pair potential. The amplitude of Re[AN(a;)] 
is suppressed due to the reduction of the proximity effect. 
On the other hand, amplitudes of Im[Ac;(a;)], Re[As(2;)], 
and Im[As(a;)] are enhanced. As regards the subdomi- 
nant components, the amplitude of Im[As(a;)] for lower 
transparent cases R = 0.75 and R— 1, is one order larger 
than those of Im[Ad(x)], and Re[As(a:)] [see Fig.lTUTb)]. 
In particular, for R = 0.75, the situation that the ampli- 
tude of Im[AN(a;)] is induced near the interface of N side 
is similar to the case oi = tt/A. Subdominant compo- 
nents which break TRS is induced also by the proximity 
effect and depends on the transparencies of the junction. 
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FIG. 9: Angle resolved LDOS at the interface iox Q = and 
R = 0. 



The corresponding LDOS for = tt/Q are plotted in 
Figs. Eta) andinib). Here, we choose Tn/T^ = 0.01, 
which is the same parameters used in Fig. 1101 The re- 
sulting line shapes of LDOS are complex reflecting on 
the complicated spatial dependence of the pair poten- 
tials. For i? = 0, we can see that the zero energy states 
are formed at the interface and the resulting LDOS has 
a ZEP due to the proximity effect. This ZEP, the origin 
of which is similar to the case of = [see Fig. Elb)], 
arises from the sign change of the pair potentials between 
Re[AN(a;)] and Re[Arf_s(a;)] through Andreev reflections 
due to proximity effect. With the increase of R, the 
LDOS has a rather broad ZEP. due to the formation of 
ABS's. This ZEP for intermediate value of R originates 
from the formation of ABS's peculiar to unconventional 
superconductors^' i.e. with sign change of the pair 
potentials between Ad(0+,O) and Ad{4>-,0). With the 
further increase of R, (see the curve with R = 0.75), since 
the imaginary components of both An {x) and Ag (x) are 
induced in the vicinity of the interface, the line shape 
of the LDOS has a complicated structure reflecting on 
the complex spatial dependence of pair potentials. For 
R — 1, the LDOS has the ZEP splitting due to the ex- 
istence of subdominant s-wave component which breaks 
TRS in the D side, i.e., Im[As(x)]. As a results, in the 
N/D junctions with 6 = tt/6, we can conclude that the 
line shapes of the LDOS around zero energy change from 
(I) the ZEP due to proximity effect, (II) the ZEP due to 
the ABS, and (III) the ZEP splitting due to the formation 
of the BTRSS, with the increase of the magnitude of R. 

Finally, we look at the angle resolved LDOS 
nY){E, (j), x) in order to understand the basic features of 
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FIG. 10: Spatial dependencies of the pair potentials for the N/D junctions with 6 = 7r/6 and Ti^/Td = 0.01. In (a), Re[AN(a;) 
Re[Ad(a;)], and in (b), Im[AN(a;)], Im[Ad(2;)], Re[As(a::)], Im[As(a;)], are plotted, respectively. 



the line shapes of iVD(£^,0). In Figs. [13 and [T3| the an- 
gle resolved LDOS n-£,{E, 4>, x) at the interface is plotted 
for 6 — n/6 and R = 0.75. In order to understand the 
role of the induced pair potential in N region, we inten- 
tionally neglect the pair potential An (a;) in the N side 
in the actual calculation of LDOS as shown in the upper 
panel of Fig. ^] and Fig. We also consider the case, 
where induced s-wave component As(x) is neglected in 
the calculation of LDOS. First we concentrate on the 
upper panels where An{x) is absent. If only the pure 
c?2.2_j^2-wave component exists in the D side, the quasi- 
particles form the ABS's at zero energy. For cj) = ±7r/3, 
since Ad((/)+, O)Ad(0_ , 0) < is satisfied, nD{E,(j),0) has 
the ZEP [see Figs. IT^ a) and ITW b)] . In the presence of 
the subdominant s-wave component which breaks TRS 
in the D side, the resulting n-[){E, (j),0) is modified. In 
this case, the line shape of njy{E, 0, 0) is not asymmetric 
around E = 0, i.e., where nDl-E, 0,0) ^ nY){—E,(j),0), 
due to the formation of the BTRSS. Therefore, the peak 
position of nD{E, cj), 0) is shifted from zero energy. On the 
other hands, for </> = ±7r/20, as Ad(0+, O)Ad(0_, 0) > 



is satisfied, riD^E, cj), 0) has gap structures [see Figs. lltjf a') 
and ll^r b)]. These results are consistent with those by 
Matsumoto and Shibai^ Next, we look at lower panels 
of Figs . [T^ and [T51 Since induced Im[AN(a;)] breaks time 
reversal symmetry, the resulting line shape of LDOS is no 
more symmetric around E — 0, even for As(a;) — 0. As 
shown in the lower panels, the presence of An (a;) forms 
a minigap around zero energy E = 0. 

In Fig. m the LDOS n^iE, (/), 0) is plotted for i? = 0, 
where the Andreev reflection only exists. In this case, the 
line shape of the LDOS can be understood only by tak- 
ing into account the pair potentials An (a;) and Ad{x). 
For <f> = 7r/3 and = ±7r/20, the LDOS has no ZEP 
since An(0+, 0)A(;((/)_|_, 0) > is satisfied. The result- 
ing nD(i?, 0, 0) has a gap structure [see Figs. EJa) and 
I14f b)]. If the quasiparticle's scattering feel different sign 
of the pair potentials between the N and D sides, we can 
expect that the nj:,{E, </>, 0) has the ZEP even for R — 0. 
For (j) = — 7r/3, since An(</>-i-, 0)Ad((/)+, 0) < is satisfied, 
then the noiE, cj), 0) has the ZEP [see Fig. lHTb)]. T hese 
situations are similar to those studied in Sec. IIII Bl 
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T^/T, = 0.01 
6=71/6 








(b) 
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e=n/6 




R = 0.75 
R= 1 



FIG. 11: The LDOS for the N/D junctions with Q = 7r/6 and 
IN/Td = 0.01. (a) i? = 0, 0.25, and 0.5, (b) R = 0.75 and 1. 



IV. CONCLUSIONS 

In this paper, we have studied quasiparticle proper- 
ties in the N/D junctions, in the presence of the prox- 
imity effect as well as the BTRSS in the D side. We 
assume the attractive interelectron potentials which in- 
duces subdominant s-wave components both in the N and 
D sides. The spatial dependencies of the pair potentials 
in the N/D junctions are determined self-consistently. 
The LDOS at the interface of the N/D junctions are 



studied in detail by changing orientational angle B and 
the transparency of the junction. Our main results are 
summarized as follows. 

i) For (110) oriented junction with Q = tt/4, the pre- 
dominant dx2-y2-waMe component in the D side is re- 
duced at the interface. The subdominant s-wave com- 
ponent which breaks TRS appears near the (110) inter- 
face in low transparent cases. In the N side, subdom- 
inant imaginary s-wave component is induced near the 
interfaces^. The resulting LDOS at the interface in the 
N side has the ZEP or its splitting depending on the 
transparency of the junction. These results are consis- 
tent with previous theoretical works liiiSiSiiSiS^ 

ii) For (100) oriented junction, the subdominant s- 
wave component becomes real number and does not 
break TRS^. With the increase of the transparency of 
the junction, the magnitude of the s-wave component in 
the N side is enhanced by the proximity effect. For fully 
transparent case, the LDOS has very sharp ZEP due to 
the formation of zero energy states. This ZEP originates 
from the fact that quasiparticles feel different sign of the 
pair potentials between the N and D sides through An- 
dreev reflection. 

iii) For 9 = 7r/6, the spatial dependencies of the 
pair potentials become very complicated. The resulting 
LDOS has a wide variety. For high barrier limit, the 
LDOS has the ZEP splitting. 

In the light of our results, although the penetration 
of the pair potential into the N side by the proximity 
effect is expected for (100) oriented junctions, the sub- 
dominant s-wave components both in the N and D sides 
arc real numbers. Thus, the BTRSS is not formed in 
the (100) junctions. The present results are inconsistent 
with previous prediction based on tunneling experiment 
on high- Tc cuprates by Kohen et al.^'^ In order to under- 
stand their experimental results, we must seek for other 
possibilities. 

In the present paper, we study proximity effect in N/D 
jmictions in the ballistic limit. In the actual junctions, we 
can not neglect impurity scattering effect. To reply this 
issue, two of the present author present theoretical works 
of charge transport in N/D junctions in the presence of 
the ABS where the N region is a diffusive metalfSLZa It 
is an interesting problem to extend the present theory in 
the diffusive regime based on the Keldysh Green's func- 
tion formalism.'''^ Although there are many works about 
ABS's in unconventional superconductor junctions up to 

71.72.73.74.75.76.77.78.79.80.81.82.83.84.85.86.87.88.89.90.91.92.93.94.95.96.97. 

proximity effects both in the presence of induced pair 
potential in the N region and the ABS's are not fully 
studied. Recent study by Lofwander is remarkable 
where the proximity effect is studied in the presence 
of both impurity scattering and induced subdominant 
component of the pair potentialsi^ Based on his detailed 
calculation, for (100) junction, a real combination of the 
proximity effect is always found. For (110) orientation, 
the s-wave component induced by the proximity effect in 
the N side breaks TRS.^-^S These results are consistent 
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FIG. 12: The angle resolved LDOS nu{E, <f>, x) at the interface with 9 = n/6 and R = 0.75. (a): LDOS without subdominant 
s-wave pair potential in the N side, (b): LDOS in the presence of subdominant s-wave pair potential in the N side. = 7r/3 
for the left panel and </> = —tt/3 for the right panel. 



with the present results. 

As a future problem, to clarify the charge transport 
property in Josephson junctions both in the presence of 

g^l04jl05jW6jl07jl08jl09jll^ 

and proximity effect is interesting. We want to extend 



the present theory towards these directions. 
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FIG. 13: The angle resolved LDOS no{E, (j), x) at the interface with 6 = 7r/6 and R = 0.75. (a): LDOS without subdominant 
s-wave pair potential in the N side, (b): LDOS in the presence of subdominant s-wave pair potential in the N side, (j) — 7r/20 
for the left panel and (p = — 7r/20 for the right panel. 
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APPENDIX A: BOGOLIUBOV-DE GENNES 
EQUATION AND QUASI-CLASSICAL 
APPROXIMATION 



Our theoretical formalism is started from the 
Bogoliubov-de Gennes (BdG) equation for unconven- 



tional spin-singlet superconductors, 

EnUnir) - Hounir) + J rfr'A(r, r')i)„(r'), (Al) 

Eninir) = -Hoinir) + J dr'A{r,r'yun{r'), (A2) 



where is the chemical potential, while t2„(r) and Vnif) 
denote the electron like and holelike components of the 
wave function 



Un{r) 
Vn(r) 

Vn{k,r) 



ikpk-r 



ikpk-r 



(A3) 



Here the quantities k and r stand for the unit vector of 
the wave number of the Cooper pair which is fixed on 
the Fermi surface (fc = fep/lfeFl), and the position of the 
center of mass of Cooper pair, respectively. After apply- 
ing the quasi-classical approximation, the BdG equation 
is reduced to the Andreev equationLSLS, 



£;„4'„(fc,r) 



A(fc,r) 



i^ivpfc- V-|-A(fc,r) 'r3 5'„(fc,r), 

(A4) 

A(fc,r) 



-A(fc,r)* 







(A5) 



The wave function ^'„(fc,r) is obtained by neglecting 
the rapidly oscillating plane-wave part following the 
quasi-classical approximationSLS. The k dependence of 
A{k,r) represents the symmetry of the pair potential. 

In the present study, we consider the case where a secu- 
larly reflecting surface or interface runs along the y direc- 
tion. Then, the pair potential depends only on x since 
the system is homogeneous along the y direction. The 
wave function ^'„(fe,r) can be written in the following 
directional notationiSiS^ 



*n(fc,a;,y) 



■\kFy\y 



(A6) 



Here ± represents the sign of the x component of the 
Fermi wave number fcpx and a(/3) = ±. We define 
a Green's function Gap{4>, x, x') and a quasi-classical 
Green's function gai3(4>,x) 



Gal3{(f>,X,x') = ^ 



ILOrii — En 



(A7) 



gap{(l),x) ±1(73)0/3 = -2h\vpx\hGal3{4',X ±0,x). 



(A8) 
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where 73 is the Pauh matrix in the directional spao 
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APPENDIX B: EVOLUTION OPERATOR 

Uai(l>a,X,x') 

The quasi-classical Green's function can be written by 
the following evolution operator Ua{4'a, x') as 

ga(3{(t),x) = Ua{(j)a,X,x')gal3{cj),x')Up'^{cj)p,X,x'), (Bl) 

where Ua((j)a, x, x') satisfies the Andreev equation 



In the above, "^^n"^ 



, x) denotes the transposition of 
^'n"'' (</'q, x). In the actual numerical calculations, we use 
UU<l^+,x,0) and t7^(0+,O,a;) given by 



U^{^+,x,0)^U^{cb+,x,0)e^^^, (B7) 
U^{cb+,0,x) - U^{4>+,0,x)e-''>. (B8) 



ih\vFx\-Q^Ua{4>a, X, x') 



-a 



iuSmh + A{(f)a,x) 



Ua{c^a,X,x'), (B2) 



with Ua{4>a,X,x) = 1. 

Hence, the evolution operators in the I side can be 
divided into a growing part and a decaying part 

f/i(0„, X, x') -AL(+)(0„, X, x')e«° 

+ AL(-'('/'c,x,x')e-''°(^-^'), (B3) 
-i.<(+)(0„,xr<(-)(0„,x')f2. 



A^) X, x') ^-^^l!"' {^.,^f^'^^'^ i^c,^') 



T2, 



where 



\VFx 



--T$^(-)(0,,a;)f2$5(+)(0„,a;) 
= const. 



(B4) 



(B5) 



(B6) 



APPENDIX C: FOR AVOIDING DIVERGENCE 
IN OUR CALCULATION 



In particular, if A'^((/)_|_, —00) [A^((/)_, 00)] ~ 0, we can 
find [C/P] 00. In this case, instead of Eq. H2U|I . we 
solve the following equation as 



d 



1 



dx \Q[,{x)^ 
+A'(,/)„,x)' 
under initial condition. 
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